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Abstract 

We consider the Newtonian limit of the theory based on the La- 
grangian 



The gravitational potential of a point mass turns out to be a combina- 
tion of Newtonian and Yukawa terms. For sixth-order gravity (p = 1) 
the coefficients are calculated explicitly. For general p one gets 



with Y^!=o c i = 1/3- Therefore, the potential is always unbounded 
near the origin. 

Wir betrachten den Newtonschen Grenzwert der durch den La- 
grangian 
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beschriebenen Theorie. Das Gravitationspotential einer Punktmasse 
hat die Gestalt von Kombinationen von Newton- und Yukawa- Termen. 
Fur Gravitationsgleichungen sechster Ordnung (p = 1) wird das Fern- 
feld von Punktmassen explizit angegeben. Fur allgemeine p wird die 
Summe der Koeffizienten vor den Yukawa- Termen durch YHI=q °i = 
1/3 beschrieben. Also ist das Potential stets unbeschrankt. 
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1 Introduction 

One-loop quantum corrections to the Einstein equation can be described by 
curvature-squared terms and lead to fourth-order gravitational field equa- 
tions; their Newtonian limit is described by a potential "Newton + one 
Yukawa term", cf. e.g. STELLE (1978) and TEYSSANDIER (1989). A 
Yukawa potential has the form exp(— r/l)/r and was originally used by 
YUKAWA (1935) to describe the meson field. 

Higher-loop quantum corrections to the Einstein equation are expected 
to contain terms of the type RU k R in the Lagrangian, which leads to a gravi- 
tational field equation of order 2£; + 4, cf. GOTTLOBER et al. (1990). Some 
preliminary results to this type of equations are already due to BUCHDAHL 
(1951). 

In the present paper we deduce the Newtonian limit following from this 
higher order field equation. The Newtonian limit of General Relativity 
Theory is the usual Newtonian theory, cf. e.g. DAUTCOURT (1964) or 
STEPHANI (1977). From the general structure of the linearized higher- 
order field equation, cf. SCHMIDT (1990), one can expect that for this 
higher-order gravity the far field of the point mass in the Newtonian limit is 
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the Newtonian potential plus a sum of different Yukawa terms. And just this 
form is that one discussed in connection with the fifth force, cf. GERBAL 
and SIROUSSE-ZIA (1989), STACEY et al. (1981) and SANDERS (1984, 
1986). 

Here we are interested in the details of this connection between higher- 
order gravity and the lengths and coefficients in the corresponding Yukawa 
terms. 



2 Lagrangian and field equation 

Let us start with the Lagrangian 

C= a k RD k R\ a p ^0. (1) 

V fc=0 / 

In our considerations we will assume that for the gravitational constant G and 
for the speed of light c it holds G = c = 1. This only means a special choice 
of units. In eq. (1), R denotes the curvature scalar, □ the D'Alembertian, 
and g the determinant of the metric. Consequently, the coefficient has the 
dimension "length to the power 2k + 2" . 

The starting point for the deduction of the field equation is the principle 
of minimal action. A necessary condition for it is the stationarity of the 
action: 



5gij 

where T y denotes the energy-momentum tensor. The explicit equations for 

5C 



are given by SCHMIDT (1990). Here we only need the linearized field equa- 
tion. It reads, cf. GOTTLOBER et al. (1990) 

pij = R ij _R g ij + 2 Y^ a k [g ij n k+1 R-U k W i i}= &nT ij , (2) 
2 k=0 



and for the trace it holds: 

77 - 1 P 

9ii ■ P ij = + 2n E a k [g ij n k+1 R] = 8nT . (3) 

n is the number of spatial dimensions; the most important application is of 
cause n — 3. From now on we put n — 3. 

3 The Newtonian limit 

The Newtonian limit is the weak-field static approximation. So we use the 
linearized field equation and insert a static metric and an energy-momentum 
tensor 

T tJ = 5°5°p, p>0 (4) 

into eq. (2). 

Without proof we mention that the metric can be brought into spatially 
conformally flat form and so we use 

9ij Vij fij i 

rjij = diag(l, -1, -1, -1) and 
fij = diag(-2$, -2*, -2*, -2*) . 

So the metric equals 

ds 2 = (1 - 2§)dt 2 - (1 + 2^)(dx 2 + dy 2 + dz 2 ) , (5) 

where $ and \I/ depend on x, y and z. 

Linearization means that the metric has only a small difference to 77^; 
quadratic expressions in and its partial derivatives are neglected. 

We especially consider the case of a point mass. In this case it holds: 
$ = $( r ), * = \J/(r), with 

r = \J ' x 1 + y 2 + z 2 , 
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because of spherical symmetry and p = m5. Using these properties, we 
deduce the field equation and discuss the existence of solutions of the above 
mentioned type. 

At first we make some helpful general considerations: The functions <3> and 
are determined by eq. (2) for % = j = and the trace of eq. (2). If these 
two equations hold, then all other component-equations are automatically 
satisfied. For the 00-equation we need R 00 : 

R 00 = -A$ . (6) 

Here is as usual 

Q2 Q2 Q2 

A = 1 1 . 

dx 2 dy 2 dz 2 

For the inverse metric we get 

g ij = diag (1/(1 - 2$), -1/(1 + 2*), -1/(1 + 2*), -1/(1 + 2*)) 

and 1/(1 — 2$) = 1 + 2$ + h($), where h(&) is quadratic in $ and vanishes 
after linearization. So we get 

In our coordinate system, / lJ equals for all i,j. For the curvature scalar 
we get 

R = 2(2A^ - A$) . (7) 

Moreover, we need expressions of the type d fc R. OR is defined by OR = 
R-ijg l \ where ";" denotes the covariant derivative. Remarks: Because of 
linearization we may replace the covariant derivative with the partial one. 
So we get 

O k R = ( l) k 2(-A fc+1 $ + 2A fc+1 ^) (8) 
and after some calculus 

- 8np = A$ + A^ . (9) 
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We use eq. (9) to eliminate \I/ from the system. So we get by help of eq. (8) 
an equation relating $ and p = m5. 

-4tt (p + 8]T a k (-l) k A k+1 p J = A$ + 6 £ a fc (-l) fe A fc+2 $ . (10) 

V fc=0 / fc=0 

In spherical coordinates it holds 

2 

A$ = -$ r + $ rr , 

r ' 

because $ depends on the radial coordinate r only. 

We apply the following lemma: In the sense of distributions it holds 

A(-e- r ' 1 ) = ^-e~ r / 1 - 4n5 . (11) 
\r / ri 2 

Now we are ready to solve the whole problem. We assume 

in I ^ \ 
$ = — 1 + ]T c * exp(-r/Z<) , Zi > . 
r \ i=0 J 

Without loss of generality we may assume U ^ lj for % ^ j. Then eq. (10) 

together with the lemma eq. (11) yield 

Svr f a k (-l) k A k ^ S = f(f + 6± a k (-l) k -^] -L-A 

i=0 V k=0 l i / 

k=0j=ki=0 H 

where t L = l\\ therefore also U ^ tj for % ^ j. This equation is equivalent to 
the system (12, 13, 14) 

X> = l/3 (12) 



i=0 



q c 



ES = 8 = 1,. ..p (13) 



i=0 i 



«r 1 + 6E«*(-l) fc *T* = i = 0,...q. (14) 

fc=0 
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From eq. (14) we see that the values U represent q + 1 different solutions of 
one polynomial. This polynominal has the degree p + 1. Therefore q < p. 



Now we use eqs. (12) and (13). They can be written in matrix form as 
( 1...1 \ ( cn \ / 1/3 \ 



1/fo • • • l/* fi 



Cl 







Here, the first q + 1 rows form a regular matrix (Vandermonde matrix). 
Therefore, we get 



j = g + l,...p 



fe=0 



with certain coefficients Xjk i.e., the remaining rows depend on the first q + 1 
ones. If Xjo ^ then the system has no solution. So A j0 = for all g + 1 < 
j < p. But for q < p we would get 



i/*? = EVi*/<? _1 
fc=i 

and this is a contradiction to the above stated regularity. Therefore p equals 
q. The polynomial in (14) may be written as 



1 l/t c 



\ 



\ i i/t p ...i/% ) 



a 



\ 



— to 



This matrix is again a Vandermonde one, i.e., there exists always a unique 
solution ( ) (which are the coefficients of the quantum corrections to 

the Einstein equation) such that the Newtonian limit of the corresponding 
gravitational field equation is a sum of Newtonian and Yukawa potential 
with prescribed lengths U. A more explicit form of the solution is given in 
the appendix. 



4 Discussion 



Let us give some special examples of the deduced formulas of the Newtonian 
limit of the theory described by the Lagrangian (1). If all the vanish we 
get of course the usual Newton theory 



$ = 



m 



A$ = -And , 



<3> and ^ refer to the metric according to eq. (5). For p=0we get for a < 



$ = 



m 



1 i p -r/V-Gai) 

+ 3 



(cf. STELLE 1978) and 



m 

q? = _ 

r 



1 p -»"/V-6ao 
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(cf. SCHMIDT 1986). For a > no Newtonian limit exists. 
For p — 1, i.e., the theory following from sixth-order gravity 
C = (R + a R 2 + ai RDRj , 



we get 



and 



where 



and 



m 

$ = — 

r 



m 

* = — 

r 



1 + c e- r / io + cie 



-r/h 



1 - c e~ r/l ° - Cl e- r/h 



1 _^ Qp 

6 2^/9a§ + 6ai 



Zo,i = \/-3a ± y9ao + 6ai . 

(This result is similar in structure but has different coefficients as in fourth- 
order gravity with included square of the Weyl tensor in the Lagrangian.) 
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The Newtonian limit for the degenerated case l = l\ can be obtained by 
a limiting procedure as follows: As we already know ao < 0, we can choose 
the length unit such that a = —1/3. The limiting case 9a 2 , + 6ai — > may 
be expressed by a\ = — 1/6 + e 2 . After linearization in e we get: 



and applying the limit e — > to the corresponding fields $ and \I/ we get 

$ = m / r {l + (1/3 + r/6)e~ r } 
= m /r{l - (1/3 + r/6)e- r } . 

For the general case p > 1, the potential is a complicated expression, but 
some properties are explicitly known, (these hold also for p — 0, 1). One gets 

$ = m/r (l + Yl Qexp(-r//i) ) 



\ i=0 / 
where J2°i — 1/3; J2 means Ya=o an d h an d c « are ( U P to permutation of 
indices) uniquely determined by the Lagrangian. 

There exist some inequalities between the coefficients Oj, which must be 
fulfilled in order to get a physically acceptable Newtonian limit. By this 
phrase we mean that besides the above conditions, additionally the fields $ 
and ^ vanish for r — > oo and that the derivatives d§/dr and d^/dr behave 
like 0(1 /r 2 ). These inequalities express essentially the fact that the U are 
real, positive, and different from each other. The last of these three conditions 
can be weakened by allowing the q to be polynomials in r instead of being 
constants, cf. the example with p — 1 calculated above. 

The equality J2 c i — 1/3 means that the gravitational potential is un- 
bounded and behaves (up to a factor 4/3) like the Newtonian potential for 



k = 1 ± ^/2tCi = 1/6 ± l/(6"\/6e) 



and 
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r PS 0. The equation $ + ^ = 2m/r enables us to rewrite the metric as 
ds 2 = (1 - 20) [(1 - 2m/r)dt 2 - (1 + 2m/r)(dx 2 + dy 2 + dz 2 )\ , 



which is the conformally transformed linearized Schwarzschild metric with 
the conformal factor 1 — 26, where 

r 

can be expressed as functional of the curvature scalar, this is the linearized 
version of the conformal transformation theorem, cf. SCHMIDT (1990). 

For an arbitrary matter configuration the gravitational potential can be 
obtained by the usual integration procedure. 



Appendix 



For general p and characteristic lengths U fulfilling < l < li < 
write the Lagrangian as 

R 



C = R 



(l 2 + ... + l 2 )R + (l 2 lf + l 2 l 2 ... + IU1 2 P )UR + 



(i 2 i(i 2 + ... + ii_jWi) nR + ■ ■ ■ + l o ■ ti ■ ■ ■ ■ • K nPR 



< l p we 



the coefficients in front of U l R in this formula read 



E 



n i i 



0<jo<ji<---<ji<P m=0 

Using this form of the Lagrangian, the gravitational potential of a point mass 
reads 



m 
r 

m 



i+i^-i^ni^-ir 1 ^ 



i=0 



7? 



i-i^-i^nil-ir 1 ^ 

i=0 jjii L i 
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For a homogeneous sphere of radius r and mass m we get 

1 p 

1 + — 12 e ~ r/h l i d i ( r o cosh(r //i) - k sinh(r /O) , 
r o i=o 

where 

e i = (-i) i+i nil-ir 1 - 

References 

BUCHDAHL, H.: 1951, Acta Math. 85, 63. 
DAUTCOURT, G.: 1964, Acta Phys. Polon. 25, 637. 

GERBAL, D., SIROUSSE-ZIA, H.: 1989, C. R. Acad. Sci. Paris, 309, Serie 
2, 353. 

GOTTLOBER, S., SCHMIDT, H.-J., STAROBINSKY, A. A.: 1990, Class. 
Quantum Grav. 7, 893. 

SANDERS, R. H.: 1984, Astron. Astrophys. 136, L21. 
SANDERS, R. H.: 1986, Astrophys. J. 154, 135. 
SCHMIDT, H.-J.: 1986, Astron. Nachr. 307, 339. 
SCHMIDT, H.-J.: 1990, Class. Quantum Grav. 7, 1023. 
STACEY, F. D. et al.: 1981, Phys. Rev. D 23, 2683. 
STELLE, K.: 1978, Gen. Rel. Grav. 9, 353. 
STEPHANI, H.: Allgemeine Relativitatstheorie, Berlin 1977. 
TEYSSANDIER, P.: 1989, Class. Quant. Grav. 6, 219. 
YUKAWA, H.: 1935, Proc. Math. Phys. Soc. Japan 17, 48. 

Remark added in proof: The cosmological consequences of these sixth- 
order field equations are recently discussed in A. Berkin and K. Maeda (Phys. 
Lett. B 245, 348; 1990) and S. Gottlober, V. Muller and H.-J. Schmidt (1990; 
preprint PRE-ZIAP 90-16). 

Received 1990 August 10 




11 



In this reprint done with the kind permission of the copyright owner we re- 
moved only obvious misprints of the original, which was published in Astronomi- 
sche Nachrichten: Astron. Nachr. 312 (1991) Nr. 2, pages 97 - 102. 

Authors's address that time: 

I. Quandt, Humboldt-Universitat zu Berlin, Unter den Linden 6, O-1086 Berlin 

H.-J. Schmidt, Zentralinstitut fur Astrophysik, 0-1591 Potsdam-Babelsberg, R.- 
Luxemburg-Str. 17a 



12 



